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Abstract 

In this work, we investigate those B c — > D ( *'T decays in perturbative QCD approach, based on 
kx factorization, where T denotes a light tensor meson. For all decays considered in this work, 
there are no contributions from factorizable emission diagrams because the emitted meson is the 
tensor meson. We find that the annihilation amplitudes are dominant in these decays due to the 
JL . large Cabibbo-Kobayashi-Maskawa elements, which are only calculable in the pQCD approach. 

The numerical results show that the predictions for the branching ratios of most decays are in the 
order of 10 -6 or even bigger, which can be observed in the ongoing experiments. We also predict 
large percentage of transverse polarizations in those W annihilation diagram dominant B c — > D*T 
decay channels. 
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I. INTRODUCTION 



After the first observation was reported in 1998 by the CDF collaboration which 
was confirmed until 2008 by CDF and DO collaboration j2| at Tevatron in excess of 5a 
significance, the study of B c meson is becoming one of the currently interesting topics, 
especially since the Large Hadron Collider (LHC) experiment ran normally. From the point 
of structure, the B c meson is a ground state of two heavy quarks' system, with a c quark 
and a b quark, which is very different from the symmetric heavy quarkonium (cc, bb) states, 
due to the flavor B = —C = ±1 carried by B c meson. Since the B c meson carries explicit 
flavor, it can not annihilate via strong interaction or electromagnetic interaction like the 
mesons consisting of cc or bb. It can only decay via weak interaction. Thus it provides us an 



ideal platform to understand the weak interaction of heavy quark flavor [3|, |4j|. Unlike the 
heavy- light B q meson (q= u, d, s), both the b and c can decay with the other as spectator, 
or they annihilate into pairs of leptons or light mesons. If more data become available, 
the B c physics must be a good place to study the perturbative and nonperturbative QCD 
dynamics, final state interactions, even the new physics beyond the standard model jjj, 4]. In 
recent years, many theoretical studies on the production and decays of B c meson have been 



done based on Operator 



roduction Expansion [5j, [6|, nonrelativistic QCD (NRQCD) and 



perturbative methods [7h 
Scora-Grinstein-Wise (ISGW) quark model 
and the perturbative QCD (PQCD) approach 



11] . QCD sum rules [1 

15-171 




131 ]. SU(3) flavor symmetry [lj 



QCD factorization approach 



18 



sgur- 



19|, 



The B meson decays involving a tensor meson have been studied in refs. 27H40I]. In 



refs. 



is 



171 ]. the authors have studied some analogous B c decays involving a tensor meson 



in final states, but only with the tensor meson as the recoiled meson. In this work, we focus 
on the B c — > D^T decays, where T denotes a light tensor meson with J p = 2 + , which is 
emitted from vacuum. We know that factorizable amplitude proportional to matrix element 
< T | j 1 * | > , where j M is the (V ± A) or (S ± P) current, does not contribute because 



28 



23, 



33 



341 ] . so these 



this matrix element vanishes from lorentz covariance considerations 
B c — > D(*>T decays are prohibited in naive factorization. To our knowledge, these decays 
are never considered in the theoretical papers due to this difficulty of factorization. In 
order to give the predictions to these decay channels, it is necessary to go beyond the naive 
factorization to calculate the nonfactorizable and annihilation diagrams. What is more, 
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the annihilation amplitudes will be dominant in considered B c — > D^T decays because 
they depend upon the large Cabibbo-Kobayashi-Maskawa (CKM) elements V c b and Ks(d)- 
It is worth of mentioning that the PQCD approach is almost the only method can do the 



quantitative calculations of the annihilation type diagrams 



41 



42|. The PQCD approach 



have successful 
B° -> D7K+ 



v 



45 



predicted the pure annihilation type decays B s — > ty + ti J43l . |44| and 



461 ] . which have been confirmed by experiments later 47|, |48[. So, for 



these annihilation dominant decays, the calculation in PQCD approach is reliable. 

In this paper, we shall study these B c — > D^T decays in the PQCD approach, which 



is based on the kj< factorization 



49M51I]. In this approach, we keep the transverse momen- 



tum of quarks, and as a result, the end-point singularity in collinear factorization can be 
avoided. On the other hand, the double logarithms will appear in QCD correction due to 
the additional energy scale introduced by the transverse momentum. Using the renormal- 
ization group equation, the double logarithms can be resumed, which results in the Sudakov 
form factor. This factor effectively suppresses the end point contribution of the distribution 
amplitude of mesons in the small transverse momentum region, which makes the calculation 
in PQCD appraoch reliable and consistent. 

In these decays, there is one more intermediate energy scale, the D meson mass. As a 
result, another expansion series of mp lmn, w ill appear. The factorization is only approved 
at the leading of mo /rriB c expansion |52|, |53J , which has also been proved by soft collinear 
effective theory (54)]. Therefore, we will take only the leading order contribution in account, 
unless explicitly mentioned. 

This paper is organized as follows. In Sec. II, we present the formalism and wave functions 
of the considered decays. Then we perform the perturbative calculations for considered decay 
channels with the PQCD approach in Sec. III. The numerical results and phenomenological 
analysis are given in Sec. IV. Finally, Sec.V contains a short summary. 



II. FORMALISM AND WAVE FUNCTION 

In order to give the predictions for these considered B c — > D^T decays, the key step is 
to calculate the transition matrix elements: 



M oc (D w T|H e// |5 c ) 



(1) 



3 



where the weak effective Hamiltonian "H e // can be written as 55] 

G 



fc q=u,c 



-V t tV tx 



10 



i=3 



(2) 



with V g 5(x) and V t b(x) (X = d, s) the CKM matrix elements. Oj (j — 1, 10) are the local 
four-quark operators: 

current-current (tree) operators 



0\ — (b a Q/3)v-A(Qi3X a )v-A, 0\ — {b a q a )v-A{<ll3Xp)y_ A , 



(3) 



QCD penguin operators 



3 = (b a X a ) v _ A J^(gJgJ) v _ A , 4 = (b a Xp) V -A ^tfptQv-A 

q> q> 

5 = (b a X a ) v _ A ^2(i? p tf p ) v+A , 6 = (l a Xp) v _ A ^2{q'pq' a 

q> q> 



V+A, 



(4) 
(5) 



electro-weak penguin operators 
3 



7 = -(b a X a ) v _ A ^2e q ,(q , pq'p) v+A , 8 = -(b a X /3 ) v _ A ^2e q >(q'^q , c 



a)V+A, 



(6) 



9 = -(6 Q X Q )v-A^ e <?'(^^)v-A, O w = -(b a Xp) v _ A ^e ql (q' p q' a ) v _ A , (7) 

q> q> 

where a and /3 are the color indices and q' are the active quarks at the scale m&, i. e. 
q' = (u,d,s,c,b). The left-handed and right-handed currents are defined as (jb a qp)v-A = 
&q7m(! - 75)<?/3 and (g^)y+A = 2a>(l + 7s) ?a respectively. The combinations a* of the 



Wilson coefficients are defined as 



56): 



ai = C 2 + Ci/3, a 2 = C x + C 2 /3, 

ai = d + a+i/3, i = 3,5,7,9, a,j = Cj + Cj_i/3, j = 4, 6, 8, 10. (8) 

In hadronic -B decays, there are several typical scales, and expansions with respect to 
the ratios of the scales are ususlly carried out. The electroweak physics higher than W 
boson mass can be calculated perturbatively. The physics between b quark mass scale and 
W boson mass scale can be included in the above Wilson coefficients Ci(/i) of the effective 



four-quark operators, which is obtained by using the renormalization group equation. The 
physics between Mb and the factorization scale is included in the calculated hard part in the 
PQCD approach. The physics below the factorization scale is nonperturbative and described 
by the hadronic wave functions of mesons, which is universal for all decay modes. Finally, 
in the PQCD approach, the decay amplitude can be factorized into the convolution of the 
the Wilson coefficients C(t), the hard scattering kernel and the light-cone wave functions 
®Mi,(B) of mesons characterized by different scales, 



A 



Tr [C{t)Q B (x u h)1>M 2 (x 2 , b 2 )$ M3 (x 3 , h)H{x u b t , t)S t ( Xi )e- s W] 



(9) 



where hi is the conjugate variable of quark's transverse momentum kir, Xi is the momentum 
fractions of valence quarks and t is the largest scale in the hard part H(xi,bi,t). The 
jet function S t (xi), which is obtained by the threshold resummation, smears the end-point 
singularities on X{ [s7|. The Sudakov form factor e~ s ^ is from the resummation of the 
double logarithms, which suppresses the soft dynamics effectively i.e. the long distance 



contributions in the large b region 



3, & Thus it makes the perturbative calculation of 



the hard part H applicable at intermediate scale, i.e., m# scale. 

In the PQCD approach, the initial and final state meson wave functions are the most 
important non-perturbative inputs. For B c meson, we only consider the contribution from 
the first Lorentz structure, like B q (q = u, d, s) meson, 



: (f + m Bc )750 Bc (x,6). 



For the distribution amplitude, we adopt the model 20]: 



■ S(x - m c /m Bc ) exp 



2i.2 
— WO 

2 



(10) 



(11) 



in which exp [— ^w 2 b 2 ] represents the kx dependence. fs c and N c = 3 are the decay constant 
of B c meson and the color number respectively. 

As discussion in ref . [27 1 , for these B c — > D^T decays, the ±2 polarizations (J p = 2 + ) 
do not contribute due to the angular momentum conservation argument. Because of the 
simplification, the wave functions for a generic tensor meson are defined by {27] 



1 

Tt 
1 



m Tf, L (pT(x) + 



T {x) + m T 



(f) s T {x) 



P ■ v 

l{x) + m T ie llupal ^e*,\n p v a (t) a T {x)\ 



(12) 



where e. = e ^ v v , and is the polarization tensor, which can be found in refs. 27-29]. The 
distribution amplitudes can be given by [27H29] 



MX 



<f) s T (x) 
<f> v T (x) 



h j. xt fr h (t) 
2V2AT 



K\x), 



2 V ^2N~ C * ± v " T1X ' J^£2N~ c dx 
The asymptotic twist- 2 and twist-3 distributions are 



(13) 



27J-|29j 



<l>\\,±(x) — 30x(l — x)(2x 



h®(x) = — (2x-l)(l-6x + 6x 2 ), h\"'(x) = 15z(l-:r)(2x-l), 



15 



(«)/ 



g [ l\x) = 20z(l-x)(2x-l), g ( l'(x) = 5(2x - l) 3 . 



(14) 



These light-cone distribution amplitudes (LCDAs) of the light tensor meson are asymmetric 
under the interchange of momentum fractions of quark and anti-quark in the SU(3) limit 



because of the Bose statistics 



28 



29]. 



For meson, in the heavy quark limit, the two-parton LCDAs can be written as refs. 



20|,|60|463] 



(D(p)\q a (z)c p (0)\0) 
(D*(p)\q a (z)c p (0)\0) 



V2NJ 



dxe lxp - z [ l5 (f + m D )0 D (x,6)] a/3 , 



V2N C 

+ MF + m D .)<i>l.(x,b)] 



x,b) 



(15) 



or the distribution amplitude for D meson, we take the same model as that used in Refs. 



61 



63] 



(j) D {x,b) 



2y/2N c 



f D 6x(l -x)[l + C D {1 - 2x)\ exp 



-co 2 b 2 



(16) 



with C D = 0.5 ± 0.1, u = 0.1 GeV and f D = 207 MeV [6j for D(D) meson and C D 



0.4 ± 0.1, u = 0.2 GeV and f Ds = 241 MeV 



64j for D S (D S ) meson. For D* meson, we take 



the same model as the D meson and determine the decay constant by using the following 
relation based on heavy quark effective theory (HQET) [651 ] . 



f D 



m 



D, 



D 



(17) 
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III. PERTURB ATIVE CALCULATION 

There are 6 types of diagrams contributing to the B c — > D^*'T decays, which are shown 
in Fig.l. The dominant factorizable emission type diagrams in most other decay modes are 
not shown here, because they do not contribute for a tensor meson emission. The second 
line are the factorizable and nonfactorizable annihilation type diagrams. 

After the perturbative calculation, the decay amplitudes for the non-factorizable emission 
diagrams in Fig. 1(a) and (b) are 

(i) (V-A)(V-A) operators: 

32 f 1 f 1 ^ 

•Menf = -^-nC F m 4 Bc I d[x] / bidbib 2 db 2 <f)B c {xi,b 1 )(f)T{x2)4 > D{x3,bi) 

Jo Jo 

X {[r D (l -X 3 ) +Xi + X 2 - 1] Eenfitajhenfix!, (1 - X 2 ), X 3 , b X , b 2 ) 

- [r D (l - x 3 ) + xi - x 2 + x 3 - 1] E enf (t b )h enf (xi,x 2 , x 3 , 61, b 2 )} , (18) 

(ii) (V-A)(V+A) operators: 

32 f 1 f 1 ^ 

■Menf = ir^C F rTm% c I d[x] / bidbib 2 db 2 (pB c {. x iibi)(t)D{.x 3 ,bi) 
<J Jo Jo 

x {[<j) s T {x 2 ){xi + x 2 + r D (x 1 + x 2 + x 3 - 2) - 1) 
+4> t T (x 2 )((xi + x 2 ){\ + r D ) - r D x 3 - 1)] 
■E enf (t a )h enf (xi, (1 - x 2 ),x 3 , 61, b 2 ) 
+ [4> t T {x 2 )(xi - x 2 + r D (xi - x 2 - x 3 + 1)) 
-(jf T {x 2 )(xi - x 2 + r D (xi - x 2 + x 3 - 1))] 

■E enf (t b )h enf (x 1 ,x 2 , x 3 , 61, 6 2 )} , (19) 
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(iii) (S-P)(S+P) operators: 

32 f 1 f 1 ^ 

•Menf = — ^-nC F m 4 Bc / d[x] / b 1 db 1 b 2 db 2 (f)B c (xi,bi)(f)T(x2)(f)D(x3,b 1 ) 
<J Jo Jo 

x [(r D (x 3 - 1) - xi - x 2 - x 3 + 2)E enf (t a )h enf (x 1 , (1 - x 2 ),x 3 , b ± , b 2 ) 

+(r D (l - x 3 ) +xi- x 2 )E enf (t b )h enf (x 1 ,x 2 , x 3 , h, b 2 )\ , (20) 

where C F = 4/3 is the group factor of SU(3) C . The hard scale t a ( b ) and the functions E en f 
and h en f can be found in Appendix A. 

Fig. 1(c) and 1(d) are the factorizable annihilation diagrams, whose contributions are 

(i) (V-A)(V-A) operators: 

Maj = 8^C F Trf B m% c dx 2 dx 3 

x {[2^ s T {x 2 )r D r T {x 3 + 1) + <p T (x 2 )x 3 ] E af (t c )h afl (x 2 , x 3 , b 2 , b 3 ) 
+ [<pT{x 2 ){2r c r D - x 2 ) + r T (-(j) t T (x 2 )(2r D (x 2 - 1) + r c ) 
+(t) s T {x 2 ){-2{x 2 + l)r D + r c ))\ E af (t d )h af2 (x 2 , x 3 , b 2 , b 3 )} , (21) 

(ii) (S-P)(S+P) operators: 



-l/A 

"3 / b 2 db 2 b 3 db 3 (f) D (x 3 ,b 3 ) 



o 



-"Mir ' ' • ■ 



Mff = -M>\ -C F f Bc nm% / dx 2 dx 3 / 4>d{x 3 M) 



JO 



x \(2(f) s T (x 2 )r T + r D (f) T (x 2 )x 3 )E af (t c )h afl (x 2 ,x 3 ,b 2 ,b 3 ) 

+ {4>T{x 2 ){2r D - r c ) + r T (<j) s T (x 2 )(x 2 - Ar D r c ) - <f> t T (x 2 )x2)) 

•E af (t d )h af2 (x 2 ,x 3 ,b 2 ,b 3 )], (22) 

with r c = m c /m Bc . m c is the mass of the c quark. t c ^, E a f and /i a /i(2) are also listed in 
Appendix A. 

The last two diagrams in Fig.l are the nonfactorizable annihilation diagrams, whose 
contributions are 

(i) (V-A)(V-A) operators: 

32 f 1 f 1 ^ 

■M-anf = — ^-C F 7rm Bc / d[x] / b 1 db 1 b 2 db 2 (f)B c (x 1 ,b 1 )(f) D (x 3 ,b 2 ) 
<J Jo Jo 

x { [<Pt(x 2 )(1 -x 1 -x 2 - r b ) - r T r D ((p t T (x 2 )(x 1 + x 2 - x 3 ) 
+<Pt( x 2)(xi + x 2 + x 3 - 2 + 4r b ))]E m f(t e )h an f 1 (xi,X2,x 3 ,b 1 ,b 2 ) 
+ [<pTi.x 2 )r D r T (-Xi + x 2 + x 3 + 4r c ) + 4 >t T {x 2 )r D r T (x l - x 2 + x 3 ) 
+<Pt(x 2 )(x 3 + r c )\ E anf (t f )h anf2 (x 1 ,x 2 ,x 3 , b u b 2 )} , (23) 



(ii) (V-A)(V+A) operators: 

32 f 1 f 1 ^ 

Manf = — T C F 7rm B c / d [ x ] / M&iM&20b c ( x i , bi)<f> D (x 3 , b 2 ) 

Jo Jo 

x { [-(4>t( x 2) + <f>T( x 2))rT(xi + x 2 - 1 - r b ) + (p T (x 2 )r D (x 3 - 1 - r b )] 
■E an f(t e )h an fi(xi,X2, x 3 , bi, b 2 ) 

+ [-(<p s T (x 2 ) + (^(^M^i - x 2 + r c ) - (f) T (x 2 )r D (x 3 - r c )] 
■E a nf(tf)h anf2 (x u x 2 , x 3 , b u b 2 )} , (24) 

with r b = m b /m Bc . i e (/), -E a «/ and h an fi( 2 ) are also listed in Appendix A. 

With the factorization formulae obtained in the above, for these B c — > DT decays, the 
total amplitudes containing the Wilson coefficients and CKM elements can be written as 

A(B C a+D°) = 9L { v: b V ud M L ^C, + v: b V cd {M L a f ai + M%jCx) 

-V t lV td [M% f (C 3 + C 9 ) + M™ f (C B + C 7 ) + M L a f{a A + a 10 ) 
+M s J(a 6 + a 8 ) + M L a ^ f (C 3 + C 9 ) + ^(C 5 + C 7 )]} , (25) 

-t#fc[A4£>(-Cs + 3a 10 /2) + MS(-C 5 + CV/2) 
+^f T f/(3C 8 /2) - M^(a 4 + aio) - M S J (a e + a 8 ) 
--^"/(Cs + C 9 ) - M L a *(C 5 + C 7 )]} , (27) 

^(S c K 2 *°L>+) = ^ {O^-M^V + -M^Cx) 

-V£V u [M% f (C 3 - C 9 /2) + MS(C 5 - C 7 /2) + M^(a 4 + a 10 ) 

(<* + ag) + M% f (C s + C 9 ) + ^(C 5 + C 7 )]} , (28) 

.A(5 C = ^^{v: b V ud M% f C 2 + v: b V cd (M%a 1 + M 

-V t * b V td [M^ f (C 3 + 2C 4 - C 9 /2 + C 10 /2) + - C 7 /2) 

+Mf n P /(2C 6 + Cs/2) + A^^(a 4 + a 10 ) + M s a f(a 6 + a 8 ) 
+-ML L /(C3 + C9) + ^(C7 5 + C 7 )]} , (29) 
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A(B C -> /•£>+) = ^ {-V t * 6 ^[^en/(^4 - C w /2) + AC/C^ - C 8 /2)]} , (30) 
.4(£ c ^a°L>+) = ^^{^y^^C 2 -^t.[A<^3(7 M /2 + A<S3a./2]},(31) 

^(S e ->^°D+) = ^ {OU^a/V + Ai^Ci) - t^Vw[A<^(a* + 010) 

+Alf;(a 6 + as) + A^S/(Cs + 0.) + ^an/(C 5 + C r )]} , (32) 

A(B c ^f!Dt) = ^l={v: b V us M^ f C 2 -V t tV ts [M^ f (2C i + C 10 /2) 

+M s e : f (2C 6 + C 8 /2))}, (33) 

A{B c ^f 2 Dt) = ^ {y c ;K s (^/V + a<S/Ci) - vsy u [M^ f (c 6 -c 7 /2) 

+M% f (C 3 + C 4 - C 9 /2 - C w /2) + M 8 * f {Ct - C s /2) 

+M^{a A + oio) + M s a f (a 6 + o 8 ) + M^ f (C 3 + C 9 ) 

+A<S/(C6 + CV)]}, (34) 

From Eq. (l4TT) . we know that 

.A(£ c -> D w / 2 ) = .A(£ c -)• £> ( *Vl) cos^ + .A(£ c -> £> w /|) sin#, (35) 
„4(5 C D^f 2 ) = A[B C ->■ DM/I) sin# - ^(5 C DW/*) cos#, (36) 

with = 7.8°. 

The amplitudes of B c — >■ D*T decay can be decomposed as 

„4(e D , e T ) = M w + i(e T D * ■ e T T *)A s + (e^n^e T D * a e^)A p , (37) 

where .4.^ contains the contribution from the longitudinal polarizations, while A s and A p 
represent the transversely polarized contributions, is the transverse polarization vector 
of D* meson, and is the vector used to construct the polarization tensors of tensor 
meson. For each decay process of B c — > D*T, the amplitudes A N , A s and A p have the same 
structures as eqs. (1251) - (134]) . respectively. The factorization formulae for the longitudinal and 
transverse polarization for the B c — > D*T decays are listed in Appendix B. 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 



The decay width of a B c meson at rest decaying into D and T meson is 



T(B C -> DT) 



•arm. 



■\A(B C ->DT)\ 



(38) 



where the momentum of the final state particle is given by 

1 



2m Br 



[ m Bc — ( m D + %) 2 ] [ m B c — ( m -D — ITIt) 2 ] ■ 



(39) 



The masses and decay constants of tensor mesons needed in the numerical calculations are 
summarized in Table [H Other parameters such as QCD scale (GeV), the mass (GeV) and 
the lifetime and decay constant of B c meson are 

= 0.25, m Bc = 6.286, f Bc = 0.489, 



r Bc = 0A6ps, u Bc = 0.6, m b = 4.8, m c = 1.5. 



(40) 



For the CKM matrix elements, here we adopt the Wolfenstein parameterization, and take 
A = 0.808, A = 0.2253, p = 0.132 and fj = 0.341 [48|. 



TABLE I: The masses and decay constants of light tensor mesons 28l. I66I. |67| 



Tensor (mass (MeV) ) 


MMeV) 


/^r(MeV) 


/ 2 (1270) 


102 ± 6 


117 ± 25 


^(1525) 


126 ± 4 


65 ± 12 


a 2 (1320) 


107 ± 6 


105 ± 21 


K|(1430) 


118 ± 5 


77 ± 14 



Like the 77 — 77' mixing, the isoscalar tensor states ^(1270) and /2(1525) also have a 
mixing and can be given by 



h = f!cos9 + flsmO, 



(41) 



with fi = ^=(uu + dd), fl = ss and the mixing angle 9 = 5.8° [68j, 7.8° |6j| or (9 ± 1)° 
48|. 
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For B c D*T decays, with three kinds of polarization amplitudes, the decay width can 
be written as 

17*1 

T(B C -> D*T) = A N | 2 +2(| A s \ 2 + | A p | 2 )). (42) 

TABLE II: Branching ratios (unit:10~ 6 ) and direct CP asymmetries (unit:%) of B c — > DT decays 
calculated in the PQCD approach . 



Decay Modes 


Class 


Br 


Adir 
A CP 


B c - 




A 


9 1 7+0.83+0.17+0.20 
z,± ' -0.71 -0.17-0.18 


a A 7+1.35 +5.33 +0.00 
'-1.15-1.59 -0.74 


B c - 


+ d°k; + 


A 


oi q+10.3 +2.81 +0.86 
0i - J -8.76 -2.86 -0.54 


n ..+0.13+0.10+0.10 
u -^-0.15 -0.22 -0.02 


B c - 


+ D+a° 2 


A 


1 ln +0.42 +0.09 -0.23 
1 - lu -0.36 -0.11 -0.26 


1 o n+4.73+10.2+0.00 
j.o.z_ 3 77 _4 .65-2 .30 


B c - 


* D+Kf 


A 


oi e+11. 3+3.10+1.01 
oi - u -9.69 -2.13 -0.63 


0.0 


B c - 


+ D+f 2 


A 


i r 1+0.58 +0.12 +0.14 
1 - dJ --0.48-0.Q9-0.16 


Q 7-1+3.45+4.09+2.70 
y -' 1 -3.97-5.21-1.59 


B c - 


+ D+f> 


A,P 


n nl 9 +0.006+0.004+0.001 
U.U1Z_ 0Q5 _ .003 -0.002 


A7 c+16.9 +10.2 +9.7 
*'- d -20.1-4.8-9.7 


B c - 


+ Df4 


C 


n nn/i7+ - 0011 +0.0016+0.0006 

U.UU<± t _ .0007 -0.0012 -0.0004 


9 n/ . +0.34 +0.62 +0.58 
z - u ^-0.37 -1.29 -0.28 


B c - 


* DfKf 


A 


1 nn+0.67 +0.20 +0.09 
x,au -0.59 -0.22 -0.07 


1 nn+°- 76 +°- 72 + 00 

i - uu -0.82 -0.50 -0.03 


B c - 


+ Dff 2 


A,P 


1 07+O.43 +0.45 +0.06 
L -° '-0.40-0.44 -0.06 


9 c-o+0.51 +1.45 +0.10 
Z - OJ -0.48 -0.72 -0.51 


B c - 


+ 


A 


AH Q+ll-9 +4.32 +1.20 
* u - y -10.7 -4.17 -0.81 


n i i +0.02+0.03+0.02 
u - iJ --0.02 -0.06 -0.00 



The CP averaging branching ratios and the direct CP asymmetries for the considered 
decay modes by using the PQCD approach are summarized in Tables [III and IIII1 The nu- 
merical results obtained from perturbative calculation are sensitive to many parameters. 
For the theoretical uncertainties in our calculations, we estimated three kinds of them: The 
first errors are caused by the hadronic parameters of mesons' wave functions, such as the 
decay constants and the shape parameters of light tensor meson, charmed meson and the 
B c meson, which are given in Sec. Eland this section. The second errors are estimated from 
the uncertainty of Aqcd — (0.25 ± 0.05) GeV and the choice of the hard scales which vary 
from 0.8t to 1.2i, which characterize the unknown next-to-leading order QCD corrections. 
The third error is from the uncertainties of the CKM matrix elements. It is easy to see 
that the most important theoretical uncertainty is caused by the non-perturbative hadronic 
parameters, which can be improved by experiments. 

It is easy to find that there are large theoretical uncertainties in any of the individual 
decay channel calculations mostly due to the shortage of the Tensor meson property. In 
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order to reduce the effects of the choice of input parameters, we define the ratios of the 
branching ratios between relevant decay modes: 

Br(B c -> D ( *1°a+) 



Br{B c ->■ D(*)+a° 2 ) 

Br(B c ->■ D^+Kf) Br(B c ->■ £>W+a!j) 



J Br(5 c ->■ £>i* )+ tf* ) 



Br(B c ^D+f 2 ) 
Br(B c ->■ D+fQ ) " 

Br(B c ^D*+f 2 ) ( 1 / /2 y cd V ! 




(43) 
(44) 

(45) 

(46) 
(47) 



Br(B c ^D*+Kf) \V2fK 5 V cs J 40' 

It is obvious that any significant deviation from the above relations will be a test of factor- 
ization or signal of new physics. 

TABLE III: Branching ratios (unit:10~ 6 ), direct CP asymmetries (unit:%) and the percentage of 
transverse polarizations i?-r(unit:%) of B c — > D*T decays calculated in the PQCD approach. 



Decay Modes 


Class 


Br 


Adir 
A CP 


R T 


B c - 


-»• D*°4 


A 


7 O/1+2.05 +0.99 +0.24 
'•^-1.75 -0.49 -0.12 


c n9 +0.54 +1.34 +0.07 
°- uz -0.54-1.37-0.51 


69.8 


B c - 


-> D*°K$ 


A 


i ci +30.1 +18.2+4.69 
±,J1 -26.5 -10.5 -3.00 


n i t-+0.02+0.05+0.03 
u - id -0.02 -0.08 -0.06 


82.5 


B c - 


-> Z?*+a° 


A 


o 7C+1.05 +0.49 +0.05 
°-'°-0.88 -0.23 -0.02 


7 QA +1.25 +4.07 +0.34 
' - y ^-l. 23-3. 87-1. 26 


68.2 


B c - 


-»• D*+Kf 


A 


1 co+30.6 +16.0 +0.00 
1JO -28.5 -14.9 -13.4 


0.0 


80.3 


B c - 


+ D*+/ 2 


A 


o qo+1.03 +0.43 +0.33 
°-°°-0.90 -0.22 -0.26 


9 /| 7+1.01 +1.55+0.82 
z -*'-l.ll-5.11-0.00 


69.7 


B c - 


+ L>*+/^ 


A 


n ncn +0.025 +0.011 +0.009 

U.U»l_ 023 _ 008 -0.009 


r fi9 +1.40 +4.63 +0.29 
d - oz -1.55 -6.30 -0.00 


45.3 


B c - 


+ 


C 


n nn^i+ - 0008 +°- 0022 + - 0006 

U.UUdl_ 00O6 -0.0015 -0.0004 


q oi +0.24 +0.52 +1.09 
°- ol -0.17 -0.81 -0.51 


12.7 


B c - 


■»• D* + Kf 


A 


q/i+1.70 +0.79+0.45 
°- y ^-1.58 -0.92 -0.28 


9 q n +0.24 +0.85 +0.01 
z - ou -0.14 -0.45 -0.01 


82.0 


B c - 




A 


cn+0.42 +0.61 +0.11 
°- uu -0.38 -0.51 -0.08 


9 nQ +0.15 +0.39 +0.10 
^• uy -0.16 -0.41 -0.40 


98.4 


B c - 




A 


inn+30.5 +19.6 +6.14 
iyu -28.1 -13.2 -3.88 


n n^fi+ - 004 +°- 011 + - 008 

U.UOD_ 003 -0.012 -0.001 


89.5 



For all considered B c — > D^T decays, the factorizable emission diagrams do not con- 
tribute, because the tensor meson can not be produced through local (V ± A) and (S ± P) 
currents. But these decays can get contributions from nonfactorizable and annihilation di- 
agrams. In fact, most of these decays are dominant by the W annihilation diagrams (A) as 
classified in the tables. There are only four decay channels, which are dominated by the color 
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suppressed (C) or penguin (P) diagrams. As we know, usually the annihilation diagrams 
are power suppressed comparing with the emission diagrams in PQCD approach. But for 
these considered decay channels, the contributions from the annihilation type diagrams are 
enhanced by the large CKM elements V cs rd) and thus play a crucial role in amplitudes. 

the predicted branching ratios are in the 



4] 



the LHC experiment, specifically the 



From Table [TT] and 1111} one can find that most o: 
order of 1CT 6 or even bigger. As stated in ref. 3, 
LHCb, can produce around 5 x 10 10 B c events each year. The B c decays with a decay rate at 
the level of 10~ 6 can be detected with a good precision at LHC experiments Jjj]. On the basis 
of our predictions, most of these B c — > D^T decays can be observed in the experiments 
soon. On the other hand, since the contributions from penguin operators are so small 
comparing with the contributions from tree operators, the direct CP asymmetries are all 
very small except B c — > D + j' 2 . For B c — > D + f 2 decay, the tree contributions from /| term are 
suppressed by the mixing angle (see 04]), to be at the same level with penguin contributions 
from /| term. The interference is sizable, thus the direct CP asymmetry is around -50%. 
Unfortunately, this decay channel is not accessible easily by current experiments due to a 
too small branching ratio. 

For B c —7- D*T decays, we also calculate the percentage of the transverse polarization 
Rt, which can be described as 



2(\A S 


2 + 


A p \ 


2 ) 


\A N 


2 + 2(\A s 


2 + 


\A p 


2 ) 



^ T I AN\2 i OH As\2 i I An\2\- (^8) 



Usually from naive factorization expectation, the longitudinal polarizations dominate the 
branching ratios of B decays. However, from numerical results shown in Table UTTl one can see 
that the transverse polarized contributions are about at the same level with the longitudinal 
polarized contributions. In fact, from eq. (IB4IIB14l) . we can find that although the transverse 
polarized contributions are power suppressed, they are also about at the same level with 
the longitudinal polarized contributions because the two factorizable annihilation diagrams 
strongly cancel with each other in the longitudinally polarized case. As a result, for these 
W annihilation diagrams dominant decays, the percentages of the transverse polarization 
are around 70% or even bigger. This large percentage can be understood as follows 71| : 
We know that the "light quark-unti-quark" pair created from hard gluon are left-handed or 
right-handed with equal opportunity. What is more, the c quark from four quark operator 



is right-handed. So the D* meson can be longitudinally polarized or transversely polarized 
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with polarization A = — 1. For the tensor meson, the anti-quark from four quark operator 
is right-handed, and the quark produced from hard gluon can be either left-handed or 
right-handed. So the tensor meson can be longitudinally polarized or transversely polarized 
with polarization A = — 1, because of the additional contribution from the orbital angular 
momentum. So the transverse polarization can become so large with additional interference 
from other diagrams. For B c — > D^ + f2, the longitudinal contributions from color suppressed 
diagrams and W annihilation diagrams strongly cancel with each other, while the transverse 
contributions can not cancel because the transverse contributions from color suppressed tree 
diagrams are too small. As a result, the ratio of transverse polarizations becomes as large 
as 98.4%. But for the color suppressed dominant B c — > D* s + a\ decay, according to the 
power counting rules in the factorization assumption, the longitudinal contributions should 
be dominant due to the quark helicity analysis 



72 



73]. The ratio is only around 10%. 



V. SUMMARY 

In this paper, we investigate B c — > D^*'T decays within the framework of perturbative 
QCD approach. We estimate and calculate the contributions of different diagrams in the 
leading order approximation of mo/mg c expansion. Most of these decays are dominant 
by the W annihilation diagrams, which are only calculable in the pQCD approach. After 
calculation, we find that the branching ratios of many decays are in the order of 10~ 6 or 
even bigger, which can be detected in the ongoing experiments. These samples of B c decays 
would provide an opportunity to study properties of B c meson and learn about the modes of 
the decays with a tensor meson emitted. Most of the direct CP asymmetries are very small 
because the penguin contributions are too small comparing with the tree contributions. We 
also predict large ratios of transverse polarizations around 70% or even bigger for those W 
annihilation dominant decays. 
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Appendix A: Related Hard Functions 



In this appendix, we summarize the functions that appear in the analytic formulas in the 
Section IHIl The first two diagrams in Fig. 1 are nonfactorizable emission diagrams, whose 
hard scales t a ^) can be determined by 



t a = max{y (xi - r 2 D )(l - x 3 )m Bc ,y\ (x 3 - 1)[(1 - rf>)(l - x 2 ) - (xi - r 2 D )\ \ m Bc , 

l/h, l/b 2 }, (Al) 



t b = max{y (xi - r 2 D )(l - x 3 )m Bc ,^\ (x 3 - 1)[(1 - r 2 D )x 2 - (x x - r 2 D )] \ m Bc , 

1/61,1/62}. (A2) 

The evolution factors E en f(t a ) and E en f(tb) in the analytic formulas (see Section IHIj) are 
given by 

E enf (t) = a s (t) exp[S Bc (t) - S T (t) - S D (t)}\ bl=h . (A3) 
The Sudakov exponents are defined as 

SbM = s (*i^f + \f ih J7<r(«.ti0). (A4) 



S D {t) = s (*3^|,&s) + 2 J* ^ lg (a s (ii)) 



(A5) 



S T (t) = s [x 2 ^^ + s^l-x^^M^j + 2 J^% q (a s (fl)), (A6) 

where the s(Q, b) can be found in the Appendix A in the ref.[50j. The function h en f can be 
given as 

h enf (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) = [9(b 2 - b 1 )K (D m Bc b 2 )I (D m Bc b 1 ) 

+8{h - b^K^DvmMUDoniBM)} 



f f [^/W\^Bb 2 ) , £> 2 <0; 
\ K (Dm B b 2 ) , ' D 2 > 0, 



(A7) 



with 



D 2 = (l-x 3 )( Xl -r 2 D ), (A8) 
D 2 = (x 3 - 1)[(1 - r 2 D )x 2 - ( Xl - r 2 D )\. (A9) 
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For the rest of diagrams, the related functions are summarized as follows: 



t c = max-{V(l - r 2 D )x 3 m Bc , l/6 2 , I/63}, 



t d = max{yx 2 x 3 (l - r 2 D )m Bc , J (I - r 2 D )x 2 + r 2 D - r 2 m Bc , l/6 2 , I/63}, (A10) 
£ a/ (t) = a s (t) ■ exp[-S T (t) - S D (t)\, (All) 



h a fi(x 2 ,x 3 ,b 2 ,b 3 ) 



Hq Jx 2 x 3 (1 - r 2 D )m Bc b 2 



6(b 2 - b 3 )H^ ] ^^F 2 m Bc b 2 j Jo L/FfmBj* 
6(b 3 - b 2 )H { 1] (^F~ 2 m Bc b 3 \ J {\[F~lm Bc b 2 



with 



h af2 (x 2 ,x 3 ,b 2 ,b 3 ) = h a fi(x2,x 3 ,b 2 ,h)\ b2 <_+ b3iF 2^ F 2, 



Ft = (1 - r 2 D )x 3 , 



(1 - r 2 D )x 2 + r 2 D - r\ 



The S t (x) is the Jet function with the expression as 57] 

2 1+2 T(3/2 + c 



S t (x) 



x(l -x)Y 



S t (x 3 ). (A12) 
(A13) 

(A14) 
(A15) 

(A16) 



v^r(i + c ) 

where c = 0.3. For the nonfactorizable diagrams, we omit the S t (x), because it provides a 
very small numerical effect to the amplitude 70 1. 



max{y £2X3(1 - r 2 D )m Bc 
1/61,1/62}, 



rj| - (1 - z 3 )(l - xi - (1 - rf,)^)^^ 



£/ = m^{Jx 2 x 3 (l-r%)m Bc ,J\r 2 + x 3 {x 1 -(l-r 2 D )x2)\m Bc ,l/bi,l/b2}, (A17) 



= a s (t) • exp[-Sfl(t) - Sr(t) - S D (t)) 



f>2=&3 ' 



(A18) 



Knfj(Xl J X 2 ,X 3 ,b 1 ,b 2 ) 



lit 
T 



0(6! - 6 2 )/7^ (Gm Bc h) Jo (Gm Bc b 2 
+6(b 2 - h)H^ (Gm B b 2 ) J (Gm Bc bi 
x , T-o 1} (\f\Gf\mBA) ,G 2 <0 



ITT TTK 

2 n 

K (Gjm Bc bi) , 



G 2 > 0, 



(A19) 
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with j = 1,2. 

G 2 = x 2 x 3 (l-r 2 D ), (A20) 

G\ = r\ - (1 - x 3 )(l - Xl -{1- r 2 D )x 2 ), (A21) 

G\ = r 2 c +x 3 ( Xl -(l-r 2 D )x 2 ). (A22) 

Appendix B: factorization formulae for B c — > D*T 

For longitudinal polarization, the decay amplitude of various diagrams and various effec- 
tive operators are 

32 /"i/A 

M enf N) = T nC F m B c / d [ X ] / M&lM&20B c ( X l > & 1 )0T (^3 , &1 ) 

<J Jo Jo 

x {[r D (l -x 3 ) -xi -x 2 + 1] E enf (t a )h enf (x 1 , (1 - x 2 ), ^3, &i, 62) 

+ [r D (l - x 3 ) + xi - x 2 + x 3 - 1] E enf (t b )h en f(xi,x 2 , x 3 , b u b 2 )} , (Bl) 

^nf N) = -^-^ c FrTm% c j d[x] j b 1 db 1 b 2 db 2 ^B c (xi,bi)(/) D (x 3 ,b 1 ) 
3 Jo Jo 

x {[4> s T {x 2 )((r D - l)(x l + x 2 ) - r D x 3 + 1) 
+0^(a; 2 )(-xi - x 2 + r D (xi + x 2 + x 3 - 2) + 1)] 
■E enf (t a )h enf (x 1 , (1 - x 2 ), x 3 , 61, 62) 
+ [-0r(x 2 )(x 2 - xi + r D (xi - X 2 - X 3 + 1)) 
+0T( a; 2)(a;2 - xi + r D (xi - x 2 + x 3 - 1))] 

■E enf (t h )h en f(x u x 2 , x 3 , b u b 2 )} , (B2) 

- A ^f^ (Ar) = ^^CW^ / / b 1 db 1 b 2 db 2 (t)B c { x iM) ( t ) T{x 2 )(l) D (x 3 ,b 1 ) 

3 Jo Jo 

x [(^d(^3 - 1) - xi - x 2 - x 3 + 2)E enf (t a )h enf (x 1 , (1 - x 2 ), x 3 , 61, 6 2 ) 

+ (r D (x 3 - 1) + xi - x 2 )E enf (t b )h enf (xi, x 2 , x 3 , bi, b 2 )\ , (B3) 

^ll(w) = g ^ -C F nf Bc m% c j dx 2 dx 3 J b 2 db 2 b 3 db 3 (f) D (x 3 ,b 3 ) 

x {[20^(x 2 )r D r T (l - x 3 ) - (j> T (x 2 )x 3 ] E af (t c )h afl (x 2 ,x 3 ,b 2 ,b 3 ) 

+ [<Pt(x 2 )x 2 + r T r c (<f> s T (x 2 ) - 0t(^2))] E af (t d )h af2 (x 2 ,x 3 ,b 2 ,b 3 )} , (B4) 
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j^sp(n) = -iQ^c F f Bc m 4 Bc -K J dx 2 dx 3 <j) D (x 3 ,b 3 ) 

x [(2<f) s T (x 2 )r T - r D (p T (x2)x 3 )E af (t c )h afl (x2,x 3 ,b2,b 3 ) 
+((<j) s T (x 2 ) - (p t T (x2))r T X2 + (pT(x 2 )r c ) 

•E a f(t d )h af2 (x 2 ,x 3 ,b 2 ,b 3 )], (B5) 

■^an\ N) = — T C F 7im B c / rf M / bidbib 2 db 2 4> Bc (xi , bi)(f) D (x 3 , b 2 ) 

Jo Jo 

x { [0t(^2)(^i + a; 2 - 1 + r 6 ) + r T r D (0^(a;2)(a;i + x 2 + x 3 - 2) 



+0t (^2 ) (zi + x 2 - x 3 ) )} E an f (t e )h an fi(x 1 ,x 2 ,x 3 ,bi,b 2 ) 

+ [-<p s T {x2)r D r T (xi - x 2 + x 3 ) + 4> t T {x2)r D r T (xi - x 2 - x 3 ) 

-<Pt(x 2 )(x 3 + r c )\ E anf (t f )h anf2 (x 1 ,x 2 , x 3 , 61, b 2 )} , (B6) 

■^ L a nf N) = --T C Frrm% c f d[x] f b 1 db 1 b 2 db 2 (t)B c { x iM) ( t ) D{x 3 ,b2) 
<J Jo Jo 

x { [-(<$rM + <^(x 2 ))r T (xi + x 2 - 1 - r b ) + (f) T (x 2 )r D (x 3 - 1 - r b )] 
■E anf (t e )h an f 1 (x 1 ,x 2 , x 3 , h, b 2 ) 

+ [-((j) s T (x 2 ) + <fr t T (x2))r T (x 1 -x 2 + r c ) - (f> T (x 2 )r D (x 3 - r c )] 
■E an f(t f )h anf 2(x u x 2 , x 3 , by, b 2 )} , (B7) 

For transverse polarization, the corresponding decay amplitudes are 

•^en/ S) = — ^Cpm^rr / d[x] f b 1 db 1 b 2 db 2 (f)B c ( x iib 1 )(j) D \(x 3 ,b 1 ) 
v3 Jo Jo 

x {[(0t(x 2 ) + 4> v T {x2)){xi + x 2 - 1)] E enf (t a )h enf (xi, (1 - x 2 ),X3,6i,6 2 ) 
+ [0t( x 2)(3;i - x 2 ) + 0y(x 2 )(-2(xi - x 2 + x 3 - l)r D + x ± - x 2 )\ 
■E enf (t b )h enf (x 1 ,x 2 , x 3 , 61, 6 2 )} , (B8) 

^C (P) = <f/ S) (B9) 

Menf S) = 7^ 7rC Fm% c / d[x] / b 1 db 1 b 2 db 2 4>B c {xi,b 1 )(j)^ ) {x 3 , bi)4%(x 2 ) 

v<j Jo Jo 

x {[r D (r D - l)(x 3 - 1)] E enf (t a )h en f(x u (1 - x 2 ),x 3 ,b 1 ,b 2 ) 

+ [rr»(rr, - l)(x 3 - 1)] E enf (t h )h enf (x u x 2 , x 3 , b u b 2 )} , (BIO) 
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M^=M^'\ (Bll) 



16 f 1 f 1 ^ 

■Menf S) = —-j=itC F m B r T J d[x] J b 1 db 1 b 2 db 2 (f>B c (xi,bi)4 £(£3 A) 

x {[<pT(x 2 ){2r D (x 1 + x 2 + x 3 - 2) - x 1 - x 2 + 1) 

+0TO r 2)O£i +a; 2 - l)]-B en /(t a )/i en /(xi, (1 - x 2 ),x 3 ,b 1 ,b 2 ) 

+ [{<Pt{ x 2) - 4>t{x 2 ))(xi - x 2 )]E enf (t b )h enf (x 1 ,x 2 ,x 3 ,b 1 ,b 2 )} , (B12) 



M%V = -M s e % s) (B13) 

M^ s) = 4\/2C F 7rf Bc r D m 4 Bc / dx 2 dx 3 / b 2 db 2 b 3 db 3 4> T D (x 3 ,b 3 ) 

Jo Jo 

x {[-M^rO^Xl ~x 3 ) + <f T (x 2 )(l + x 3 ))] E af (t c )h afl (x 2 , x 3 ,b 2 ,b 3 ) 
+ [r T ((l>T(x 2 )(x 2 - 1) + <^(x 2 )Or 2 + 1)) 

~(j) T (x 2 )r c ] E af (t d )h af2 (x 2 ,x 3 ,b 2 ,b 3 )} , (B14) 



Kf v) = M L J is) (B15) 

Mff s) = 8>/2C F 7r/B c m| c / dx 2 dx 3 / b 2 db 2 b 3 db 3 <j) T D (x 3 ,b 3 ) 

Jo Jo 
x {M^tO^) + 0t(^2))] E af (t c )h afl (x 2 , x 3 ,b 2 ,b 3 ) 

- [r D ((t)l(x 2 )(rl - 1) + 2<f) v T (x 2 )r T r c )] E af (t d )h af2 (x 2 ,x 3 ,b 2 ,b 3 )} , (B16) 

M S ^ P) = M s a ^ s) (B17) 

16 Z" 1 f 1 ^ 

M tnf S) = -^ C F^r D m Bc J d[x] J b 1 db l b 2 db 2 (f) Bc (x 1 ,b 1 )(f)^(x 3 ,b 2 ) 

x { [-<p 1 r(x 2 )r D (x 3 - 1) - 20^(x 2 )r T r b ] E anf (t e )h anfl (x 1 , x 2 , x 3 , b u b 2 ) 

+ [<pT(x2)rDX 3 + 2(p v T (x 2 )r T r c ] E anf (t f )h anf2 (x 1 ,x 2 ,x 3 ,b 1 ,b 2 )} , (B18) 



M L a ^ p) = M L a L J s) \^r T , (B19) 
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M 



LR(s) 
anf 



16 „ „ f 1 , r , /- 1/A 



C F 7rm| c / d[x] / b 1 db 1 b 2 db 2 (f)B c (xi,b 1 )(f)^(x 3 ,b2) 



Jo Jo 



x {[-(^(^2) + (p v A x 2))r T (xi +x 2 -l-r b ) 



+(pT(x 2 )r D (x 3 - l-r b )] E anf (t e )h anfl (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) 



M^tM + <t> v T (x 2 ))(xi -x 2 + r c ) 



+(j>T(x 2 )r D (x 3 - r c )] E anf (t f )h anf2 (x 1 ,x 2 ,x 3 ,b 1 ,b 2 )} , 



(B20) 




(B21) 
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